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Subject:  An  Analysis  to  Facilitate  the  Interpretation  of  Boundary 

Layer  Measurements 

Abstract:  A  computer  aided  boundary  layer  analysis,  which  draws  on 

available  theoretical  and  empirical  data,  is  presented. 

The  analysis  allows  boundary  layer  data  to  be  extracted 
from  measurements  influenced  by  a  normal  pressure  gradient, 
comparison  of  boundary  layer  measurements  against  standard 
laminar  and  turbulent  profiles,  and  evaluation  of  important 
boundary  layer  parameters.  This  technique  is  illustrated 
by  example. 
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Nomenclature 

Law  of  the  wall  constant  (=  5.0) 

Skin  friction  coefficient 

Error  between  the  data  and  the  wall-wake  equation 

Similarity  function 

Clauser's  shape  factor 

First  shape  factor 

Second  shape  factor 

Reference  length 

Number  of  data  points 

Number  of  data  points  in  the  inviscid  region 
Maximum  number  of  data  points  that  could  possibly 
be  in  the  inviscid  region 
Static  pressure 

Static  pressure  at  the  boundary  layer  edge 

Reference  Reynolds  number 

Momentum  thickness  Reynolds  number 

Streamwise  velocity 

Boundary  layer  velocity 

Inviscid  velocity 

Measured  composite  velocity 

Dimensionless  velocity  in  the  inner  boundary  layer 
Shear  or  friction  velocity 
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Velocity  at  the  boundary  layer  edge 
Reference  velocity 
Normal  velocity 

Coles'  universal  wake  function 
Streamwise  coordinate 

Coordinate  normal  to  the  blade  surface 
Dimensionless  coordinate  normal  to  the  blade  surface 
in  inner  boundary  layer  variables 

Falkner-Skan  streamwise  pressure  gradient  parameter 

Clauser's  equilibrium  parameter 

Boundary  layer  thickness 

Energy  thickness 

Displacement  thickness 

Defect  thickness 

Dimensionless  normal  similarity  variable 
Momentum  thickness 

Von  Karman's  mixing  length  parameter  (=  0.41) 
Kinematic  viscosity 

Dimensionless  streamwise  similarity  variable 
Coles'  wake  parameter 
Fluid  density 
Wall  shear  stress 


Streamfunction 
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Nomenclature  [continuation] 

Subscripts 

i  The  data  point 

wall  At  the  wall 

C  Derivative  with  respect  to  C 

Superscripts 

'  Derivative  with  respect  to  n 
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Introduction 

Over  the  past  several  years,  detailed  boundary  layer  velocity  profiles 
have  been  made  using  laser  Doppler  velocimetry  in  the  periodic  two-dimensional 
flow  over  a  double  circular  arc  compressor  blade  in  cascade.  The  initial 
incidence  angle  was  5”  so  that  laminar,  transitional,  turbulent  and  separated 
velocity  profiles  were  all  encountered.  Two  typical  profiles,  one  on  the 
suction  surface  and  one  on  the  pressure  surface,  are  shown  in  Figures  1  and  2. 
The  effect  of  a  strong  normal  pressure  gradient  is  apparent  and  even  the  most 
fundamental  properties  of  the  boundary  layers  are  masked.  Moreover, 
experience  has  shown  that  even  without  the  complication  of  a  strong  normal 
pressure  gradient,  the  state  of  a  boundary  layer  (e.g.,  laminar  or  turbulent) 
is  often  not  self  evident  and  the  method  for  calculating  the  properties  of 
the  layer  are  often  not  obvious. 

In  the  following,  a  method  of  boundary  layer  analysis,  which  draws  on 
available  theoretical  and  empirical  formulations,  is  presented.  The  analysis 
details  methods  of  reconstructing  boundary  layer  profiles  from  data 
influenced  by  a  normal  pressure  gradient,  of  comparing  data  against  standard 
laminar  and  turbulent  profiles  (with  or  without  strearawise  pressure 
gradients),  and  of  generating  important  boundary  layer  parameters  (e.g., 
thicknesses  or  shear  velocity).  The  analysis  is  complete  for  laminar, 
transitional,  and  turbulent  boundary  layers.  Extension  of  the  techniques 
to  separated  profiles  and  wakes  is  underway. 

Normal  Pressure  Gradient 

Surface  curvature,  blowing,  or  suction  may  induce  significant  streamline 
curvature  in  a  flow.  Streamline  curvature  in  turn  causes  a  normal  pressure 
gradient  which  results  in  a  cross  streamline  gradient  in  the  inviscid  velocity 


-6- 


25  May  1985 
WCZ:SD:lhz 


profile.  As  in  Figures  1  and  2  then,  the  freestream  velocity  does  not  reach 
a  constant  value  (the  edge  velocity,  Boundary  layer  analyses,  however, 

assume  that  the  location  of  the  boundary  layer  edge  is  known  and  that  a 
constant  freestream  velocity  exists  outside  the  boundary  layer.  To  use  or 
formulate  a  boundary  layer  analysis,  a  procedure  must  be  used  to  account  for 
the  effect  of  the  normal  pressure  gradient. 

One  method  which  accounts  for  these  effects  was  outlined  by  Mellor  and 
Wood  [1971]  and  Ball,  Reid  and  Schmidt  [1983].  The  method  assumes  that  the 
measured  velocity  profiles  represent  composite  profiles.  This  implies  that 
each  of  the  profiles  has  a  region  where  the  viscous  effects  predominate,  a 
region  in  which  viscous  effects  are  negligible,  and  an  intermediate  region 
in  which  the  viscid-inviscid  results  match.  Mathematically,  the  measured 
composite  profile  is  the  sura  of  a  boundary  layer  profile  and  an  inviscid 
profile,  less  what  appears  in  both.  The  last  quantity  is  commonly  called 
the  edge  velocity,  Ug.  That  is 

u  =u,.„+u,  -U  . 
meas  b£  inv  e 

Clearly,  both  the  boundary  layer  velocity  u^ji  and  the  measured  velocity  u^jeag 
must  go  to  zero  at  the  wall,  so  that 

U  =  (u.  ) 

e  inv^ygll 

and  the  scheme  reduces  to  finding  the  value  of  the  edge  velocity  at  the  wall. 

The  validity  of  applying  this  procedure  to  velocity  data  is  difficult  to 
establish.  However,  it  does  allow  for  a  consistent  method  to  analyze  measured 
velocity  profiles  provided  that  the  inviscid  region  can  be  properly  identified 
from  the  data.  There  does  not  appear  to  be  a  rigorous  way  to  do  this,  so  that 
as  in  the  prior  reported  study  (Deutsch  and  Zierke  [1984]),  a  consistent 
method  which  produces  plausible  results  is  adopted. 


-7- 


25  May  1985 
WCZ : SD : Ihz 


A  least  squares  technique  is  used  to  fit  a  polynomial  to  the  inviscid 
velocity  profile.  Since  the  inviscid  profiles  may  have  significant  curvature, 
the  option  of  fitting  a  linear,  quadratic,  or  cubic  polynomial  is  provided  in 
the  analysis.  The  difficulty  lies  in  choosing  the  data  points  to  be  included 
in  the  least  squares  analysis.  To  minimize  this  problem,  the  number  of  data 
points  used  in  the  polynomial  fit  was  varied.  First,  a  maximum  number  of 
points,  N  ,  which  could  possibly  be  within  the  inviscid  region  was 

tI13X 

determined.  For  3p/ 3y  >  0  (see  Figure  1),  N  is  taken  to  be  the  number 

max 

of  points  between  the  point  furthest  from  the  wall  and  the  point  of  maximum 

velocity.  For  3p/ 3y  <  0  (see  Figure  2),  the  data  point  at  which  the  profile 

slope  changes  by  at  least  50Z  is  used  instead  of  the  point  of  maximum 

velocity.  For  the  data  in  Figure  1,  N  was  determined  to  be  19,  while 

^max  found  to  be  24  for  the  data  of  Figure  2.  Many  applications  have 

shown  that  for  the  region  of  0.55  N  <  N.  <  0.95  N  ,  the  value  of  U 

max  inv  max  e 

is  relatively  independent  of  the  number  of  points  in  the  fit.  This  is 

clearly  shown  for  the  data  of  Figures  1  and  2  in  Table  1.  Each  polynomial 

was  extrapolated  to  the  wall  to  obtain  and  a  mean  and  a  standard  deviation 

of  all  values  of  U  were  determined.  The  degree  of  the  polynomial  was  chosen 
6 

to  minimize  the  standard  deviation  of  U  .  These  standard  deviations  have 

e 

been  observed  to  be  quite  small  (-  0.5%).  A  u  profile  was  calculated 

D 

using  the  mean  value  of  U  and  a  polynomial  fit  of  u.  .  For  u 

e  "  inv  inv’ 

N.  =  0.75  N  was  chosen  as  the  number  of  data  points  in  the  fit  as  the 
inv  max 

found  in  this  way  was  consistently  close  to  the  average  U^.  A  smoothed 
spline  fit  of  the  boundary  layer  velocity  profile  was  used  to  calculate  the 
boundary  layer  thickness,  5.  5  is  taken  as  the  position  at  which 
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u,  0  =  0.99  U 
bJl  e 

In  Figures  3  and  4,  the  reconstructed  boundary  layer  profiles  (triangles) 
corresponding  to  Figures  1  and  2  are  shown.  As  might  be  anticipated,  the 
effect  of  the  normal  pressure  gradient  does  not  penetrate  far  into  the 
boundary  layer. 

Deutsch  and  Zierke  [1984]  discuss  the  plausibility  of  this  approach  by 
considering  a  comparison  of  the  shape  of  the  resulting  turbulence  intensity 
versus  normalized  distance  (y/5)  from  the  wall  plots  against  classical 
measurements.  They  conclude  that  the  technique  does  indeed  give  reasonable 
results.  It  should  be  noted  that  an  alternate  method  of  accounting  for  the 
effects  of  the  normal  pressure  gradient  using  the  equations  for  the  boundary 
layer  parameters  was  developed  by  Kiock  [1983].  That  technique  was  found 
here,  however,  to  be  too  difficult  Co  employ. 

Preliminary  Boundary  Layer  Analysis 

Comparisons  between  boundary  layers  are  often  made  in  terms  of  simple 
integral  thicknesses  or  their  ratios  (shape  factors). 

In  order  to  calculate  the  boundary  layer  integral  parameters,  the 
velocity  data  points  must  be  fit  with  a  mathematical  curve.  A  parametric 
cubic  spline  was  used  for  the  curve  fit.  This  curve  parametrically  develops 
the  u^^  and  y  ordered  pairs  as  independent  functions  of  the  overall  arc 
length  of  the  curve.  Even  this  curve  fit  tends  to  oscillate  when  there  is 
a  cusp  in  the  profile.  Therefore,  a  smoothing  routine  was  added  to  the 
spline  fit.  Near  Che  wall,  a  parametric  cubic  spline  was  fit  between  the 
data  point  nearest  to  the  wall  and  the  zero  velocity  point  that  would  occur 
at  the  wall.  For  boundary  layer  profiles  not  well  resolved  in  the  near  wall 
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region,  this  latter  fit  represents  the  largest  potential  error  in  the 
calculation  of  the  integral  parameters.  Shown  in  Figures  5  and  6  are  the 
calculated  spline  fits  to  the  boundary  layer  profiles  of  Figures  3  and  4. 


Boundary  Layer  Integral  Parameters 

The  most  common  of  the  parameters  which  characterize  the  boundary  layer 
are  integral  thicknesses  based  on  the  conservation  of  mass,  momentum  and 
energy.  In  effect,  the  boundary  layer  acts  to  displace  the  streamlines  in 
the  flow  outside  the  boundary  layer  away  from  the  wall.  A  displacement 
thickness,  d*,  can  be  defined  as  the  distance  by  which  the  solid  surface 
would  have  to  be  displaced  to  maintain  the  same  mass  flowrate  in  a 
hypothetical  inviscid  flow.  For  steady,  incompressible  flow,  6*,  can  be 
defined  as 


6*  -  /  (1  -^)dy  . 

0  e 

<5  is  probably  the  most  fundamental  integral  thickness  and  is  commonly  used 
as  a  normalizing  factor  in  presenting  data  since  the  alternative  the 
boundary  layer  thickness,  is  difficult  to  measure.  A  momentum  thickness, 

9,  can  be  determined  from  the  steady,  incompressible  momentum  equation  as 


9= 


hi 


u 

0  e 


Jdy 


The  momentum  thickness  represents  the  momentum  loss  due  to  the  presence  of 
the  boundary  layer  and  is  proportional  to  the  drag  when  no  streamwise 
pressure  gradient  exists.  The  momentum  thickness  is  used  in  many  empirical 
correlations.  Finally,  a  steady,  incompressible  energy  thickness  can  be 


defined  as 
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r  ^ 

U  ^  „2 

0  e  U 


jdy 


Reynolds  numbers  can  be  formed  based  on  all  three  of  these  integral 
thicknesses  using  the  boundary  layer  edge  velocity.  They  are  especially 
useful  in  describing  the  boundary  layers  just  before,  during,  and  just 
after  transition. 

Some  parameters  describe  the  shape  of  the  velocity  profile.  The  first 
and  second  shape  factors  of  the  velocity  profile  are  defined  as 

”l2  =  f- 

and 


Note  that  the  definitions  require  that  both  shape  factors  be  greater  than 
unity.  For  laminar  boundary  layers,  the  first  shape  factor  lies  between 
3.5  and  2.3.  Transition  brings  about  a  considerable  drop  in  giving  a 
turbulent  boundary  layer  in  which  lies  between  1.3  and  2.2.  Laminar 
separation  takes  place  at  a  value  of  near  3.5  while  turbulent  separation 
takes  place  at  a  value  of  near  2.2. 

In  che  current  study,  the  integral  parameters  and  shape  factors  were 
found  by  integrating  the  spline  fit  using  a  trapezoidal  rule  with  very  fine 
spacing.  Values  of  these  parameters  and  the  associated  Reynolds  numbers 
are  given  in  Table  2  for  each  of  the  profiles  shown  in  Figures  5  and  6. 

Based  on  the  work  of  Purtell,  Klebanoff  and  Buckley  [1981]  and  Murlis,  Tsai, 
and  Bradshaw  [1982],  one  would  anticipate  that  the  Re^  for  the  profile  of 
Figure  6  is  too  low  to  support  turbulence. 
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Similarity  Solutions  for  Laminar  Boundary  Layers 

For  a  constant  streamwise  pressure  gradient,  laminar  boundary  layers  can 
be  plotted  in  such  a  way  as  to  collapse  them  all  to  the  same  curve.  Such 
laminar  boundary  layers  are  said  to  be  self-similar  and  solutions  to  the 
boundary  layer  momentum  equation  may  be  solved  using  a  similarity  solution. 

A  similarity  solution  reduces  the  number  of  variables  in  the  equation  by  using 
a  coordinate  transformation.  The  boundary  layer  momentum  equation  for  steady 
flow  can  be  writteii  as 


u  ^  +  IH.  =  _  i.  ^  9^u 

3x  3y  p  dx  -  2  * 

3y 

Several  different  coordinate  transformations  have  been  used  for  similarity 
solutions.  One  such  transformation  is  the  Levy-Lees  transformation  (see 
Cebeci  and  Smith  [1974])  which  is  given  as 


and 


where 


K  = 


x/£ 

/  (■ 


0 


U 

e 


ref 


1 


U. 


)(|) 


£U 


Re  = 


ref 


Again,  these  transformation  variables  are  for  an  incompressible  flow  with  a 

constant  value  of  viscosity.  The  variables  U  ^  and  I  refer  to  a  reference 

ref 

velocity  and  a  reference  length,  respectively.  A  similarity  function,  f, 
can  then  be  defined  from  the  strearafunction,  ij;,  as 
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where 


and 


V(x,y)  =  f(?,n) 


The  boundary  layer  momentum  equation  for  steady,  incompressible  flow  can  now 
be  transformed  into  the  following  differential  equation 


where 


f'"  +  ff"  +  e[i  -  f’^]  =  25(f’f’ 


6 


U  dC 
e 


S  is  a  streamwise  pressure  gradient  parameter  and  must  be  determined  from  the 
measured  streamwise  static  pressure  distribution.  These  measurements  are  also 
required  to  determine  the  values  of  When  f  has  a  pria  i or  a  superscript, 
it  refers  to  a  derivative  with  respect  to  n.  A  derivative  with  respect  to  C 
is  noted  by  an  C  for  a  subscript. 

A  similarity  solution  which  is  quite  useful  was  developed  by  Falkner 
and  Skan  [1931].  Their  solution  is  for  laminar  boundary  layers  flowing  over 
wedges  where  the  streamwise  pressure  gradients  are  constant.  In  terms  of 
the  Levy-Lees  coordinates,  the  Falkner-Skan  solution  involves  no  changes 
in  the  ^-direction.  Thus  all  derivatives  with  respect  to  ?  are  zero  and 
the  boundary  layer  momentum  equation  for  steady,  incompressible  flow  is 
reduced  to  the  following  ordinary  differential  equation 


f 


t  I  t 


+  ff "  +  6[l  -  f =  0  . 
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The  solution  of  this  equation  is  found  in  the  present  boundary  layer  analysis 
from  the  solution  scheme  of  Hoffman  [1983].  The  Falkner-Skan  solution 
represents  a  very  good  approximation  of  the  laminar  boundary  layer  solution 
in  two  cases:  first,  when  the  streamwise  pressure  gradient  changes  only 
slowly  in  the  streamwise  direction,  and  second,  near  the  leading  edge  of  the 
boundary  layer  surface  (since  5  is  very  small  and  the  Levy-Lees  transformed 
equation  is  approximated  quite  well  by  the  Falkner-Skan  equation).  In  the 
equation  above  when  B  is  equal  to  zero,  the  Falkner-Skan  solution  is  reduced 
to  the  original  similarity  solution  of  Blasius  [1908]  for  laminar  boundary 
layers  with  no  streamwise  pressure  gradient.  Separation  of  the  laminar 
boundary  layer  will  occur  for  a  S  with  a  value  of  -  0.199, 

The  profile  of  Figure  6  is  shown  compared  against  a  Falkner-Skan  profile 
at  an  equivalent  6  in  Figure  7.  Although  the  streamwise  pressure  gradient  is 
not  constant  here,  the  comparison  is  quite  good.  The  profile  of  Figure  6 
(via  Figures  2  and  4)  is  laminar.  Integral  parameters  and  skin  friction 
values  can  be  computed  from  the  Falkner-Skan  solution.  These  are  presented 
in  Table  2;  where  appropriate  they  are  compared  to  the  values  calculated  from 
the  spline  fit. 

Wall-Wake  Velocity  Profile  for  Turbulent  Boundary  Layers 

Turbulent  boundary  layers  are  commonly  divided  into  different  regions. 
The  innermost  region  is  dominated  by  viscous  shear  and  is  self-similar  for 
all  turbulent  boundary  layers.  This  region  is  called  the  viscous  sublayer 
and  is  described  by 


u 

T 


V 
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or 

+  + 
u  =  y 

where  u"*"  and  y"*"  are  called  inner  variables,  is  called  the  shear  or 

friction  velocity.  Outside  of  the  sublayer  but  still  very  close  to  the  wall, 
the  velocity  is  logarithmic  with  distance  as 

=  7  l-n(y'^)  +  C  . 

The  viscous  sublayer  and  logarithmic  region  overlap  in  a  region  called  the 
buffer  layer.  Collectively  the  sublayer,  buffer  layer  and  logarithmic  layer 
are  called  the  "law  of  the  wall."  Streamwise  pressure  gradients  have 
essentially  no  effect  on  this  region.  Outside  of  this  wall  region,  the 
streamwise  pressure  gradients  are  important  and  the  velocity  profile  exhibits 
a  wake-like  form.  Coles  [1956]  developed  an  equation  for  the  wake  region 
called  the  "law  of  the  wake."  His  composite  equation  included  a  wake-like 
function  added  to  the  logarithmic  equation.  This  wall-wake  equation  can  be 
written  as 

=  7  lti(y'^]  +  C  +  W(^) 

where 

W(-^]  =  2  sin^(||)  =  1  -  cos  (-2^)  . 

W  is  called  Coles'  universal  wake  function  and  is  normalized  to  be  zero  at 
the  wall  and  has  a  value  of  W  =  2.0  at  y  =  6.  Coles'  wake  parameter  11  brings 
in  the  effect  of  the  streamwise  pressure  gradient  and  has  a  value  of 
approximately  0.5  for  a  zero  pressure  gradient  flow.  Coles  and  Hirst  [1968] 

later  examined  the  data  from  a  number  of  turbulent  boundary  layer  experiments 

and  determined  that  the  von  Karman  mixing  length  parameter,  <,  should  be 

0.41  and  the  law  of  the  wall  constant,  C,  should  be  5.0. 
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The  analysis  of  turbulent  boundary  layers  should  include  a  fit  of  the 
data  to  the  wall-wake  equation.  A  least  squares  fit  of  the  data  to  the 
equation  was  chosen  with  u^  and  IT  as  the  variables  to  be  determined  and  6 
considered  to  be  a  known  quantity.  The  error  between  each  data  point  and 
the  wall-wake  equation  is 


u  y .  u  u  II 

^  [1  - 


Now  the  minimum  error  squared  of  all  the  data  points  must  be  found  with 
respect  to  both  u^  and  H.  Taking  the  partial  derivative  of  the  error 
squared  for  each  data  point  with  respect  to  u^  and  H  yields 


db.  ,  y.u  ,  „  Try. 

=  2E,  {-  ln(-^)  +1+0+  -  [l-  cost-^)]} 


u  Try 

^  =  2E.{^  [l  -  cos(-^)]}  . 


Summing  for  all  the  data  points  and  setting  the  two  expressions  equal  to  zero 


gives  two  equations 


N  9E^ 

V  _ L 

^  9u 
i=l 


N  as? 

^  lrr  =  o 
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which  can  be  solved  simultaneously  for  u^  and  11  to  give  the  minimum  squared 
error.  A  secant  method  is  used  to  solve  the  simultaneous  non-linear 
equations.  A  similar  method  was  presented  by  Sun  and  Childs  [1976]  for 
compressible  turbulent  boundary  layers  using  the  velocity  transformation 
of  van  Driest  [1951]. 

In  order  to  use  the  fit  of  the  wall-wake  equation,  the  range  of  data 
points  to  be  used  in  the  fit  must  be  considered.  White  [1974]  states  that 
the  logarithmic  region  does  not  hold  for  y^  <  35  (corresponding  roughly  to 
y/6  <  0.02).  If  the  wake  component  is  large.  Coles  and  Hirst  [1968]  state 
that  the  region  being  fit  should  not  include  y/ 6  >  0.9.  They  suggest  that 
this  number  should  be  reduced  to  0.75  for  a  zero  streamwise  pressure  gradient 
and  0.6  for  a  vanishing  wake  component.  Here  the  region  0.05  <  y/6  <  0.75 
was  used;  sensitivity  to  the  specific  range  of  data  points  was  not  large. 

The  data  of  Figures  5  (via  Figures  1  and  3)  is  shown  in  Figure  8  in 
inner  variables.  Parameters  determined  from  the  wall-wake  fit  are  given  in 
Table  2. 

Streamwise  pressure  gradients  can  have  a  strong  effect  on  the  outer 
region  of  a  turbulent  boundary  layer.  Typical  profiles  at  an  Reg  of  5,000 
for  values  of  II  of  -  0.5,  0.5,  2.0  and  5.0  are  shown  in  Figure  9(a)  for 
illustration.  We  should  note  that  the  value  of  Re„  controls  the  size  of 
the  lograithmic  region.  Shown  in  Figure  9(b)  are  the  computer  profiles  for 
Reg  of  500,  1,000,  5,000  and  10,000  and  a  n  of  0.5.  Note  the  difficulty  in 
discerning  a  logarithmic  region  for  RCg  of  500. 

Most  analyses,  short  of  direct  computation,  consider  boundary  layers  in 
which  the  pressure  gradient  is  constant,  so  that  the  profiles  are  self-similar 

with  downstream  distance  —  so  called  equilibrium  layers.  Clauser  [1954,1955] 
conceived  of  a  parameter  as 
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6 


c 


x* 

0  e 

T  dx 
w 


to  characterize  equilibrium.  Here  is  the  wall  shear  stress.  A  turbulent 
boundary  layer  with  a  constant  8^  has  outer  region  similarity  and  is  called 
an  equilibrium  turbulent  boundary  layer.  All  of  the  gross  properties  of  that 
boundary  layer  can  be  described  with  a  single  parameter.  Clauser  [1954,1955] 
chose  the  parameter  G,  where 


<=  -  i  * 

0  T 


and 


U  -  u,  . 

4  ■  /  (-  u— -jdy  • 

0  T 


G  is  called  Clauser 's  shape  factor  and  A  is  termed  the  defect  thickness. 
Values  of  G  and  A  for  the  profile  of  Figure  5  are 


and 


G  =  19.47 

A  =  142.1  ram 


The  cascade  profiles  are  non-equilibrium.  In  practice,  one  expects  non- 
equilbrium  to  be  the  rule  rather  than  the  exception. 


Skin  Friction  Coefficient  of  a  Turbulent  Boundary  Layer 

An  empirical  expression  for  the  skin  friction  coefficient  of  a 
turbulent  boundary  layer  was  developed  by  Ludweig  and  Tillman  [1949].  The 
curve-fit  expression  from  their  experimental  data  is 
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—0.268  —0.678  Hit 

C.  -  0.246  Re^  10 
t  o 


Skin  friction  coefficient  (c^j,  wall  shear  stress  and  friction  velocity 

(u^)  can  be  related  by 


T 


The  friction  velocity  calculated  from  Ludweig  and  Tillman's  expression 
agrees  well  with  that  calculated  from  the  wall-wake  fit  for  all  the  non- 
separated  turbulent  profiles  measured. 


Transition 

Transition  represents  the  region  in  which  a  laminar  boundary  layer 
becomes  turbulent.  The  length  of  the  transition  region  depends  strongly 
on,  among  others,  the  streamwise  pressure  gradient  and  the  freestream 
turbulence  Intensity.  The  incomplete  transition  of  the  pressure  surface 
profiles  is  shown  in  Figure  10.  Here  the  transition  is  termed  incomplete 
as  the  profile  at  97.8%  chord  does  not  contain  a  logarithmic  region;  also, 
Reg  is  below  that  believed  capable  of  sustaining  full  turbulence.  Note 


the  comparison  with  the  Falkner-Skan  profiles;  particularly  the  thickening 
of  the  profile  in  the  near  wall  region. 
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Conclusions 

A  computer  aided  analysis  of  boundary  layer  data  is  described.  The 
analysis  allows  for 

(1)  The  reconstruction  of  boundary  layer  profiles  from  data  influenced  by 
normal  pressure  gradients  (both  and  6  are  determined), 

(2)  A  calculation  of  boundary  layer  integral  parameters  through  a  smoothed 
spline  fit  of  the  data, 

(3)  (a)  A  comparison  of  the  data  against  a  laminar  Falkner-Skan  profile 
at  equivalent  pressure  gradient  parameter  (B)  and/or 

(b)  A  curve  fit  to  the  law  of  the  wall/wake  of  Coles  (u_^  and  II  are 
determined),  and 

(4)  A  calculation  of  skin  friction  from  Ludweig-Tillman' s  equation. 

The  analysis  has  been  successfully  used  Co  interpret  both  cascade 

boundary  layers  (air)  and  the  boundary  layer  on  an  axisymmetric  microbubble 
body  (water).  An  extension  to  separated  profiles  and  wake  profiles  is 


planned. 
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Table  1.  Variation  of  Ug  and  6  with  the  Number  of  Parts  used  in  the 
Polynomial  Fit  of  the  Inviscid  Region. 


Suction  ! 

Surface  Boundary  Layer  at  53.6% 

Chord  (3p/3y  >  0) 

^inv 

Ug,  m/sec 

6,  cm 

11 

33.64 

1.130 

12 

33.74 

1.142 

13 

33.83 

1.154 

14 

34.00 

1.177 

15 

33.84 

1.155 

16 

33.76 

1.146 

17 

33.64 

1.132 

33.78  ±  0.13 

1.148  ±  0.016 

Pressure 

Surface  Boundary  Layer  at  57.2% 

Chord  (3p/3y  <  0) 

^inv 

Ug,  m/sec 

6,  cm 

14 

23.37 

0.231 

15 

23.38 

0.233 

16 

23.39 

0.234 

17 

23.41 

0.237 

18 

23.40 

0.235 

19 

23.40 

0.235 

20 

23.40 

0.235 

21 

23.40 

0.235 

22 

23.40 

0.235 

23.39  ±  0.02 

0.234  ±  0.002 
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Table  2.  Boundary  Layer  Parameters. 


Parameters 

Suction  Surface  Boundary 
Layer  at  53.6%  Chord 
(3p/3v  >  0) 

Pressure  Surface  Boundary 
Layer  at  57.2%  Chord 
(3p/3y  <  0) 

Spline  Fit 

Fit  of  Wall- 
Wake  Eq . 

Spline  Fit 

Falkner-Skan 

Solution 

5*,  (cm) 

0.350 

0.346 

0.083 

0.059 

6  ,  (cm) 

0.180 

0.182 

0.029 

0.022 

(cm) 

0.288 

0.291 

0.046 

— 

Re^* 

7881 

7795 

1295 

835 

CD 

4041 

4095 

451 

310 

6489 

6557 

213 

— 

«12 

1.95 

1.90 

2.87 

2.70 

«32 

1.61 

1.60 

1.58 

— 

u^,  (m/sec) 

0.850^^^ 

0.823 

0.385^^^ 

0.337 

n 

/.  /.O^ 

A 

7 

c 

B 

— 

— 

— 

-  0.0588 

(1)  Calculated  from  the  Ludweig-Tillman  Equation. 
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Figure  3.  Reconstructed  boundary  layer  velocity  profile  (9p/3y  >  0) . 


Falkner-Skan  Approximation 


Figure  8.  VJall/wake  approximation  to  the  boundary  layer  velocity  pro 
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